Dynamical instabilities in density-dependent hadronic relativistic models 
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Unstable modes in asymmetric nuclear matter (ANM) at subsaturation densities are studied in the 
framework of relativistic mean-field density-dependent hadron models. The size of the instabilities 
that drive the system are calculated and a comparison with results obtained within the non-linear 
Walecka model is presented. The distillation and anti-distillation effects are discussed. 
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I. INTRODUCTION 

Many efforts are recently being done in order to under- 
stand the supernova evolution. Of particular interest is 
the scenario in the aftermath of a core bounce, where a 
large number of neutrinos is produced and radiated out 
towards the infalling matter from the outer layers onto 
the core. The mean free path of the neutrinos and their 
interaction with matter can be an explanation for the 
mantle ejection during the explosion. In [l[ the opacity 
of the nuclear non-uniform neutron-rich matter is cal- 
culated in a semiclassical approach to describe neutrino 
scattering. In the present work we will investigate some 
general properties of the non-uniform matter at the crust 
of a compact star within different relativistic models. In 
particular we will investigate the dynamical collective un- 
stable modes and study the isospin content of the non- 
homogeneous phase of asymmetric nuclear matter. 

In two previous works 0, 0] we have investigated the 
influence of the electromagnetic interaction and the pres- 
ence of electrons on the unstable modes of npe matter at 
zero and finite temperature within the NL3 parametriza- 
tion of the non- linear Walecka model (NLWM) A\. This 
parametrization describes the ground-state properties of 
both stable and unstable nuclei. 

Models with density-dependent meson-nucleon cou- 
plings are an alternative approach for the description 
of nuclear matter and finite nuclei Non- linear self- 
interactions of the mesons in constant coupling models 
are substituted by density-dependent meson-nucleon cou- 
pling parameters and are motivated by Dirac-Brueckner 
calculations of nuclear matter. 

The parametrization introduced by Typel and Wolter, 
which we will refer as TW [(|, describes finite nuclei 
and nuclear matter with similar quality as non-linear 
parametrizations and has a more reasonable extrapola- 
tion to extreme conditions: high density and large charge 
asymmetry. In [7j a parametrization denoted DD-ME1 
used the same density dependence of TW for the coupling 
parameters, but adjusted the parameters in a different 
way. More recently the parametrization DD-ME2 Q has 
been developed as an improvement of DD-ME1 in order 
to obtain better fittings to excitation energies of isoscalar 



monopole and isovector dipole giant resonances. Other 
possibilities for density dependent parametrizations are 
found in the literature 

In [10 ] it was shown that the thermodynamical insta- 
bilities at subsaturation densities of NLWM parametriza- 
tions with constant couplings differ from the behavior 
of relativistic nuclear models with density-dependent pa- 
rameters. In particular, in the last models the distillation 
effect is not so strong and follows more closely the be- 
havior of non-relativistic nuclear models. In the present 
work we will study the dynamical instabilities within 
density-dependent relativistic models (DDRM) and will 
compare them with the results obtained with the NL3 
parametrization of NLWM. 

This investigation will be performed in the framework 
of the Vlasov formalism p], [H, 0- We will study the 
role of isospin and the modification of the distillation 
phenomenon due to the presence of the Coulomb field 
and electrons. 

In section II we review the Vlasov equation formalism 
for nuclear neutral matter including electrons and the 
electromagnetic field. In section III the dispersion rela- 
tion is displayed and in section IV the numerical results 
are shown and discussed. Finally, in the last section the 
most important conclusions are drawn. 



II. THE VLASOV EQUATION FORMALISM 

We start from the lagrangian density of the relativistic 
TW model |6| including electrons interacting with the 
electromagnetic field 
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where fi M „ = d^Vv — dyV^, = d^A v - d v A^ and 
= d^K - d^bf, - r p (6 M x b v ). 
The parameters of the model are: the nucleon mass M, 
the electron mass m e , the masses of the mesons m s , m v , 
nip, the electromagnetic coupling constant e = ^/47r/137 
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and the density-dependent coupling parameters T s , T v 
and r p , which are adjusted in order to reproduce some 
of the nuclear matter bulk properties, using the following 
parametrization: 



r i(p) = Ti{paat)9i(x), i = s,v 



with 



9i(x) 



l + b i (x + dj)' 2 
l + c t (x + di) 2 



where x = pj p sa t and 



T p (p) = T p (p sat ) exp[-a p (x - 1) 



(2) 



(3) 



(4) 



In the sequel we will present results obtained with TW 
and DD-ME2. The values of the parameters rrii, Ti, m, 
bi, Ci and di, i — s,v,p are given in Table |U 



NL3 [4] TW [6] DD-ME2 [8] 



B/A (MeV) 


16.3 


16.3 


16.14 


Po (fm- 3 ) 


0.148 


0.153 


0.152 


K (MeV) 


272 


240 


250.89 


£ sym . (MeV) 


37.4 


32.0 


32.3 


M'/M 


0.60 


0.56 


0.572 



TABLE II: Nuclear matter properties. 



and the corresponding one-body hamiltonian h 
diag(/i p , h n , h e ) , with 

hi = ^(v-V l f + M* 2 + V 0l , i = p,n 

and 



TW [6] DD-ME2 [8] 



m s (MeV) 


550 


550.1238 


m v (MeV) 


783 


783.0000 


m p (MeV) 


763 


763.0000 


r s (p sa t) 


10.72854 


10.5396 


Tu (Psat ) 


13.29015 


13.0189 


r p (p sa t) 


7.32196 


7.3672 


a s 


1.365469 


1.3881 


b s 


0.226061 


1.0943 


Cs 


0.409704 


1.7057 


d 3 


0.901995 


0.4421 


a v 


1.402488 


1.3892 


b v 


0.172577 


0.9240 


Cy 


0.344293 


1.4620 


d v 


0.983955 


0.4775 


a p 


0.515 


0.5647 



TABLE I: Parameters of the density-dependent models. 



Notice that in these density-dependent models the non- 
linear terms are not present, in contrast with the usual 
non-linear Walecka model (NLWM). For comparison we 
summarize in Table HI the nuclear matter properties at 
saturation calculated for the models we will use. For 
the NL3 parametrization of the NLWM the lagrangian 
density has the same structure as (JTJ) plus the non-linear 
terms, namely 

£ = £{gs,gv:9p) - ^k</> 3 - ^A^ 4 , 

where the meson- nucleon coupling constants g s , g v , g p re- 
place r s ,r„,r p and k, A are the self-coupling constants 
for the non-linear terms. 

In order to determine the time evolution of the sys- 
tem we introduce the one-body phase-space distribu- 
tion function in isospin space /(r, p, t) = diag (f p , f n , f e ) 



h e = v(p + eA) 2 + m 2 - eA , 

where M* = M — T s cf) denotes the effective baryon mass 
and 



> 



Vo* - T V V + ^ nbo + eAo^p- + S fl 



V i = V v V + ^T i b + eA—^ 



with Ti = 1 (— 1) for protons (neutrons). The last expres- 
sions contain the contribution of a rearrangement term 
given by 



Op^^^-^ 



due to the density dependence of the coupling parameters 

The time evolution of the distribution function is de- 
scribed by the Vlasov equation 



^ + {fi,hi} = 0, 



P, n, e, 



(5) 



where {, ) denotes the Poisson brackets. It has been ar- 
gued in [3, EH that ([5|) expresses the conservation of 
the number of particles in phase space and is, therefore, 
covariant. 

The equations of motion for the fields are obtained 
from the Lagrangian and are given by 



d 2 4> 

dt 2 



V 2 + m20 = r s p s (r,t) 



(6) 



d 2 V» 

dt 2 

d 2 b^ 



2 V 



r v f(v,t)+d»(d u vn , (7) 



nrl V¥+m^ = ^(M) + ^(V), (8) 
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— -V 2 A» = e[#(r,i)-#(r,t)] , (9) 
where the scalar density is 

Ps (r,t) = 2j2 f^M^Pit) — - (10) 



2— p,n 



Si 



The components of the baryonic four-current density 
are 



j (r,t) = 2^ JAPffapj) =p P + Pn , (11 



i— p,n 



(12) 



i—p,n 



where p p , p n are the proton and neutron densities. The 
electron four-current density has components 



JO,:' 



r >*)= 2 E / (03/e(r,P,O=Pe 



J.('.*) = 2E /(2^ /e(r ' P ' 







p + eA 



(13) 



(14) 



where p e is the density of electrons and the components 
of the isovector four-current density are 
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,o(r,t) = 2 £ 



d 3 p 
(2^)3 



/i(r,p,i)Ti = p p - p n , (15) 



and is described by the distribution function /o(r, p) = 
diag(9(P 2 p -p 2 ), 6(P 2 „-p 2 ), 6(P 2 e -p 2 )) and by 
the constant mesonic fields (defined with a (0) super- 
script) which obey the following equations m 2 4>^ = 
r o (0) m 2 V {0) - T 7 (0) V (0) - m 2 6 (0) - ^7 (0) 

& f = 0, A® = 0, and = 0. 

Collective modes in the present approach correspond to 
small oscillations around the equilibrium state, and they 
are described by the linearized equations of motion [111 ]. 
We take for the distribution function f = fo + Sf and, 
as in [ill ] we introduce a generating function S(r, p, t) — 
diag (S p , S n , S e ) , defined in isospin space such that the 
variation of the distribution function is 



Sfi - {SiJoi} = -{S h p 2 }S(P Fi -p 2 ) 



(17) 



In terms of this generating function, the linearized Vlasov 
equations for Sfi are equivalent to the following time evo- 
lution equations 



^ -I- {S ei h Qe } = 5h e = ~e 
at 



dt 



+ {Si, h ai \ = Shi 



SA - 

M* 
> 



p-<5A 



£0e 



SV oi 



(18) 



P-SV 1 



(19) 

i = p, n, where we have taken linear variations for the 
fields. In equation (jT8|l So e = \Jp 2 + m 2 and in equation 



OH hoi = vV +M* 2 + V 0l = Eo + Voi . The linearized 
equations of the fields are obtained using the procedure 
already presented in [l3| . 



III. SOLUTIONS FOR THE NORMAL MODES 
AND DISPERSION RELATION 



j 3 (r,*)=2^ /^ /l(r ' P ' i)2 ^ Tl ' (16) 



with Ei = \J (p — Vi) 2 + M* 2 , i — p,n and e e = 

y/Jp + eA) 2 + m 2 . 

The four-currents j^, jg and j% satisfy the continuity 
equations |l3j c^j M = 0, d^j^ = and d^j^ = 0. 
Substituting these continuity equations into {2} and (jHJ, 
the following relations between the components of the 
vector mesonic fields are obtained: 

rofoV" = fd„T v , m 2 p d^ = jgd^Tp. 

For constant coupling parameters the above relations re- 
duce to the usual relations between the components of a 
vector field 

d^V = 0, d^V = 0. 

At zero temperature, the state which minimizes the 
energy of asymmetric nuclear matter is characterized 
by the Fermi momenta Pfi, * = P,n, Pfc = Pfp 



The longitudinal normal modes of the system, with 
momentum k and frequency to are well described by the 
ansatz 

Ti = ^i.^exp [i(wt - k ■ r)] 

for the fields and 

Sj(r, p, t) = 5i(cos6»)exp [i(wt - k • r)] , j = e, p, n, 

for the generating functions, where 9 is the angle between 
p and k. A different choice of the generating function 
would allow to study the transverse modes [lj]. This, 
however, will not be carried out in the present work. For 
the longitudinal modes 6V* = SV* = , Sb* = 6b% = 
and SAl = SAl = . 

Equations (p~8j) and (fT9|) are written in terms of the 
amplitudes A^i related to the transition densities by 
Spi = fp^rPoiAoi, and they read 

1 + FPPL p FP n L p C p A e L p \ ( A^p \ 
F n P L n 1 + F nn L n A un = 0, 

C e /L e \A^ e ) 

(20) 
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with A ui = x S U i(x) dx, Li = 2 — s t In ffj^j^ where 



Si = Lu/uJoi = u)/(kVFi), — tjt" being the Fermi 
velocity of particle i, Efi z 



M* 



p,n, e Fe = 



\J Pp e + m e 2 . We also have 



F ij = \pV Ws _Qij Wp _ GvWv 
+ ^— J {GsDi + TiGpD + GvD 



M* dT, 



bo dr r 



£ Fj Op 2 Op 

+ \{1 + ^{1 + ^)0% 



p2 



-^e Fl , (21) 



and 



with 



, 2 1 pf 



Fj 



2tt 2 fc 2 V Fl 



2tt 2 (w 2 



j = s, p, v, and ui 2 = fc 2 
ft 2 + m 2 , with fv? .* = 



L s,effi v 



, = k 2 +m 2 vl uj 2 = 



T 2 S (dp s /dM*) . All the 



other quantities are defined in the Appendix. 

From (|2"D|) we get the following dispersion relation 

[1 - Cf L e ] [1 + L P F PP + L n F nn + L p L n (F pp F nn 



_p P n pnpy _ C e lC v lL e L v {\ + L n F nn ) = 0. (22) 

In order to study the instabilities of the system, we look 
for solutions of the dispersion relation with imaginary 
frequencies. These modes are obtained by replacing s 
with i[3 in the expression for Li. 



IV. RESULTS AND DISCUSSION 

In the present section we compare the dynamical spin- 
odals, direction of instability and most unstable modes 
obtained with NL3, TW and DD-ME2. For reference, 
in Fig. [T]we compare the symmetry energy of the three 
models and the /3-equilibrium equation of state (EoS) 
at low densities. It is known that NL3 symmetry en- 
ergy grows nearly linearly with density and has a quite 
high value at saturation in comparison with TW and DD- 
ME2. Thus in /3-equilibrium matter the proton fraction 
increases very quickly and reaches values which allow for 
the direct URCA process, and therefore predict a too fast 
cooling of neutron stars, already at density values close 
to the saturation density. The two models we will con- 
sider with density-dependent couplings have very similar 
symmetry energies and predict similar proton fractions. 



While npe matter is thermodynamically stable within 
TW and DD-ME2 models, for NL3 there is still a small 
unstable region 0. Even being thermodynamically sta- 
ble, npe matter is unstable with respect to perturbations 
with certain wavelengths. The region of instability is lim- 
ited by the spinodal surface which, for a given k transfer, 
is obtained from the dispersion relation and corresponds 
to the surface on which the eigenmode is zero. In Fig. [2] 
we plot the spinodal for three values of k: 11, 75 and 150 
MeV. The value k = 75 MeV defines, except for small cor- 
rections, the envelope of the spinodals for k values. In Fig 
[2K) we include the spinodals obtained in three different 
situations: only neutron-proton (np) matter excluding 
the Coulomb field felt by the protons, together with np 
matter and neutron-proton-electron (npe) matter includ- 
ing the Coumlob interactions. While the first situation 
is not realistic, but allows a comparison with the ther- 
modynamical limit, the second describes neutron-proton 
matter and the third one, stellar matter. For fc = 75 
and 150 MeV the results for np matter with Coulomb in- 
teraction essentially coincide with those for npe. Again 
for np matter, now with no Coulomb field for fc = 11 
MeV, the results practically coincide with the thermody- 
namical spinodal (which corresponds to fc = MeV). At 
fc = 150 MeV the effect of the electrons and the Coulomb 
field is very small, as expected from the 1/fc 2 behavior of 
the Coulomb field 0. 

The following conclusions may be taken: for fc = 11 
MeV the spinodal for np matter with Coulomb is much 
smaller than the corresponding spinodals for npe matter 
due to the attractive force between protons and electrons 
in the last case; for npe symmetric matter (p p = p n ) the 
three models considered have similar results but differ- 
ences occur for asymmetric matter, DD-ME2 showing 
instabilities at larger densities for the largest asymme- 
tries. 

In FigO)) we include the /3-equilibrium EoS for npe 
neutrino- free matter Y v — and for npe^ matter as in su- 
pernovae with a constant lepton fraction Yj, = Y e + Y v = 
0.4 [l^|. The crossing of these EoS with the spinodal tell 
us that there is a non-homogeneous region in the star, 
at low densities. The density at the inner edge of the 
crust, as predicted by the present calculation, is given 
in Table IIIII For neutrino trapped matter the values 
shown are only an upper limit because for T ^ the 
size of the instability region is smaller. In this situation, 
the three models give similar results because the matter 
considered has a very high proton fraction (y p ~ 0.3), 
therefore closer to symmetric matter, where parameter 
sets are expected to coincide. However, for neutrino free 
matter, the density value at the inner edge of the crust is 
very sensitive to the model because we are dealing with 
highly asymmetric matter where the largest differences 
between models arise. 

We next analyze the direction of the instability defined 
by the ratio of the fluctuations, Sp p /Sp n , corresponding 
to the eigenmode that becomes imaginary. In Fig [3^) 
and b) we plot 5p p /Sp n as a function of fc for two proton 
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FIG. 2: Dynamical spinodals for a) k = 11 MeV, b) k = 75 MeV, c) fc = 150 MeV 



model Y V = QY L = 0.4 
NL3 0.050 0.082 
TW 0.076 0.084 
DD-ME2 0.073 0.083 

TABLE III: Density at the inner edge of the crust of a com- 
pact star 



fractions y p = 0.1, typical of neutrino free stellar mat- 
ter, and yp — 0.35 which, as quoted above, would be 
found in stellar matter with trapped neutrinos, and for 
two densities, p = 0.15 po and 0.3 po. We include, for ref- 
erence, a dashed thin line which indicates the correspond- 
ing Pp/pn ratio. In Fig. [3J;) we fix k and for the same 
proton fractions referred above we show the dependence 
of Spp/Sp n on the density. Some conclusions are in order: 
at low densities, the distillation effect, which corresponds 
to Sp p /Sp n > pp/p ni is similar for all the models, with 
NL3 slightly less efficient for larger asymmetries. This 
can also be observed from Fig [3};) and was also seen in 
the thermodynamical instability calculations at low den- 
sities [10|,|l7|. For larger densities, both DDRM are less 
effective than NL3 in the reposition of symmetry. In- 
deed, from Fig [3k) we clearly observe that for densities 
larger than the ones considered in Fig. [3^) and b) NL3 
becomes the model which more efficiently describes the 
distillation effect, while TW and DD-ME2 keep showing 



a behavior which is similar among themselves, and almost 
independent of the density. The differences between the 
two types of models are larger for larger asymmetries. 

For the np calculation with no Coulomb field the 
Spp/Sp n ratio is almost constant with respect to k, 
though slightly increasing, specially for small values of 
y p . In addition to it, if the Coulomb field is included this 
ratio becomes much smaller than for the no Coulomb 
case, crossing even the p p /p n line for k < 25 MeV, for 
the largest proton fraction considered here. This is the 
anti-distillation effect already discussed in 0. In Table 
IIVI we show, for several pairs of asymmetry-density, the 
maximum k values for which the anti-distillation occurs 
in the three models. All models have similar values al- 
though they are slightly larger for NL3. These values 
are never very large: we get k < 25 MeV and decreasing 
values of k with increasing proton fraction. As discussed 
before, the distillation effect is larger for matter with no 
electrons, for in this situation protons do not couple to 
the electrons. 

This effect may have important consequences in stel- 
lar matter. In fact, in a supernovae explosion 99% of the 
energy is carried away by the neutrinos. Neutrinos inter- 
act strongly with neutrons (large weak vector charge of 
the neutron) and therefore the way neutrons clusterize 
is important to determine the neutrino mean free path. 
Neutrinos may couple strongly to the neutron-rich matter 
low-energy modes present in this explosive environment 
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50 100 150 200 250 50 100 150 200 250 02 04 06 08 1 

k(MeV) k(MeV) ' p\fm~ 3 ) 



FIG. 3: Direction of instability as a function of the momentum transfer for p = 0.15 po and 0.3 po for (a) y p = 0.1 and (b) 
y p — 0.35 and as a function of density (c) for k = 100 MeV and y p = 0.1, 0.35, for np matter only. 




0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 0.02 0.04 0.06 0.08 0.1 

P (far 3 ) p (far 3 ) p (far 3 ) 



FIG. 4: Most unstable modes: a) growth rates and b) associated size of clusters for asymmetries y p = 0.1, 0.3, 0.5 and the 
models NL3, TW and DD-ME2. Results including electrons (thick lines) are compared with np matter results with no Coulomb 
interaction (thin lines). 



and revive the stalled supernovae shock. 

The system is driven to the non-homogeneous phase by 
the mode with a larger growth-rate. In Fig. 2]we plot the 
growth-rate of the most unstable mode as a function of 
density for np matter without Coulomb interaction and 
npe matter. The wavelength associated with these modes 
is related to the size of the inhomogeneities formed. In 
Fig. HJd) half of the wavelength, which corresponds to 
the size of the clusters formed, is plotted as a function 
of density for the proton fractions and models considered 
in Fig. 5k). As expected, in all the cases the presence 
of electrons reduces the growth-rate and the size of the 
clusters; this effect is more pronounced for larger densi- 



ties. 

For very small densities (p < 0.1 po) all three models 
have a similar behavior, characterized by a large growth- 
rate. As density increases, all the curves have similar 
slopes, but considerable differences between the mod- 
els arise. For symmetric matter TW behaves like NL3 
with the largest values for the growth rate and the size 
of the associated clusters. As asymmetry increases TW 
still maintains the largest instability, but NL3 changes its 
behavior becoming closer to DD-ME2 with the smallest 
growth-rate. 

The size of the instabilities that drive the system is of 
the order of 4 — 10 fm. For small k the unstable mode 
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Vp P 


NL3 TW DD-ME2 


0.1 0.15/90 


6 




0.3p 






0.25 0.15p 


15 13 


12.7 


0.3p 


14 13.5 




0.35 0.15po 


25 22.4 


21.8 


0.3p 


24.8 24.5 


24.6 



TABLE IV: Maximum value of k in MeV for which the anti- 
distillation effect occurs. 



disappears due to the quenching of the instability: p- di- 
vergence of the Coulomb energy. In the large k limit the 
effect of the Coulomb interaction goes to zero. Larger dif- 
ferences between NL3, TW and DD-ME2 occur at densi- 
ties and proton fraction of interest for /3-equilibrium stel- 
lar matter. In particular TW predicts larger clusters at 
densities above ~ 0.02 fm~ 3 . The size of the clusters cal- 
culated agree with the results of a density functional with 
relativistic mean-fields coupled with the electric field [l8[ . 

V. CONCLUSIONS 

We have investigated the low densities instabilities in 
density-dependent relativistic hadronic models (DDRM) 
and compared them with previous results obtained 
within NLWM, namely with the NL3 parametrization. 

The spinodal region shows that both DDRM used here 
present instability regions larger than NL3, except for 
large k in which case DD-ME2 has a smaller spinodal 
region. These differences occur mainly at larger isospin 
asymmetry. From the astrophysical point of view, this 
could mean differences in low density stellar matter, 
namely the crust properties of compact stars. In par- 
ticular, we have seen that while the predicted inner crust 
edge density for stellar matter with trapped neutrinos 
is very similar in all models considered, in cold stellar 
matter with no neutrinos the differences are large. For 
DDRM this density is about 50% larger than the one for 
NL3. 

It was shown that except for the lowest values of den- 
sity, density-dependent parametrizations predict lower 
distillation effects. At low densities this trend is no more 
true with DDRM showing results which are similar to 
NL3 or even slightly larger for small proton fractions. 

For small k an anti-distillation effect is present in npe 
matter and np matter with Coulomb interaction. It is 
for the NL3 parametrization that this behavior sets on 
at larger k values and it is present in this model even for 
very large asymmetries (see Table IIV[) . This will have 
an important effect on the scattering of neutrinos which 
escape the proto-neutron star: a large neutron fraction 
implies a larger weak force interaction. 

We have predicted the formation of clusters with sizes 



ranging from ^ 4 fm to 10 fm. These limits depend on 
the proton fraction and larger clusters are formed in more 
asymmetric matter. 

We finally conclude that different parametrizations of 
DDRM have similar properties and different from other 
models with constant couplings. Their predicitve power 
will depend on their ability of satisfying constraints both 
having astrophysical origin or laboratory measurements 
0- 

Neutrino opacity plays a crucial role but it is not 
the only mechanism of energy accounting in a super- 
nova. The plasmon-decay into neutrino-antincutrino 
pairs should also be considered in neutron star evolu- 
tion. In [2(J, we have studied plasmons in stellar matter 
within constant-coupling relativistic models. Plasmons 
are currently being studied in DDRM, as we have shown 
in [2l| . There, nuclear plasmon modes were found at zero 
temperature. We are also carrying out finite temperature 
calculations and we expect that this contribution will al- 
low estimations of neutrino production, due to neutrino- 
antineutrino decay, in nuclear matter under neutron star 
conditions. 
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APPENDIX: Dispersion relation coefficients 

The expressions used in Eq. [21] read as: 



P 



= r„ 



M* 



( dp s 



e Fi - \dM* J Q \ dp J 

dr. 



G p \i — 

G p 2i = 

G v 
G v i 



(o) . 

d P Jo' 
Gpij G p 2i, 

r ■ dT 



nT p [ l 



G V \G V 2, 



= T„ 



,(0) 



9 P Jo 
' k 2 



dp 



m 2 ,; 



9 P 7o 



(0) 
P3 i 



mi 



G V 2 — 



GsDi 
G p D 
G V D 



= H, 



dps 



b 2 i d£s 

*\dM*J \dpJ > 



An, I \k) Tp { d'p)/^' 



dT, 



7i?, v k ) v 



dT v 
dp 



M* 



+ 



V 



b 



dT s 
d P /o 

(0) 



-,(°) 



,(o) 



d 2 T t 
dp 2 



dr L 



■4 0) 



d 2 r v 

d P 2 /o 

5 2 r 



dp j ™ V 3p 2 



In the previous expressions the zero in the superscripts 
and subscripts on p and derivatives, respectively, mean 
that they are calculated with respect to the static back- 
ground on which the oscillations take place. 
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